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ABSTRACT 

Using new Hubble Space Telescope Faint Object Spectrograph, New Technology 
Telescope EMMI, and Keck HIRES spectra of the gravitationally-lensed double QSO 
HE 1104-1805 AB ( 

z em — 2.31), and assuming UV photoionization by a metagalactic 
radiation field, we derive physical conditions (ionization levels, metal abundances and 
cloud sizes along the lines of sight) in five ClV+Mgil absorption systems clustered 
around z = 1.66 along the two lines of sight. Three of these systems are associated 
with a damped Lya (DLA) system with logAT(Hl) = 20.85, which is observed in the 
ultraviolet spectra of the bright QSO image (A). The other two systems are associated 
with a Lyman-limit system with log AT (Hi) = 17.57, seen in the fainter image (B). The 
Civ and Mgn line profiles in A resemble those in the B spectra, and span Av ~ 360 
km s . The angular separation 9 = 3.195" between A and B corresponds to a 
transverse proper separation of S± = 8.3 h^Q kpc, for qo = 0.5 and a lens at z = 1. 

Assuming that the relative metal abundances in these absorption systems are the 
same as observed in the DLA system, we find that the observed iV(Clv)/iV(Mgii) 
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ratios imply ionization parameters of logT = —2.95 to —2.35. Consequently, these 
clouds should be small (0.5-1.6 kpc with a hydrogen density kh <^ 0.01 cm -3 ) and 
relatively highly ionized. The absorption systems to B are found to have a metallicity 
0.63 times lower than the metallicity of the gas giving rise to the DLA system, 
Zdla ^ 1/10 Z & . 

We detect Ovi at z = 1.66253 in both QSO spectra, but no associated Nv. Our 
model calculations lead us to conclude that the C IV clouds should be surrounded 
by large (~ 100 kpc) highly ionized low-density clouds (ran ~ 10~ 4 cm" 3 ), in which 
O vi, but only weak Civ absorption occurs. In this state, logT > —1.2 reproduces the 
observed ratio of N(0 Vl)/iV(N v) > 60. 

These results are discussed in view of the disk/halo and hierarchical structure 
formation models. 



Subject headings: galaxies: abundances — gravitational lensing — quasars: absorption 
lines — quasars: individual (HE 1104-1805 AB) 



1. Introduction 



The physical nature of the ionized gas observed in high redshift QSO absorption systems 
has not been clearly determined so far. One possible mechanism capable of ionizing this gas is 
ultraviolet (UV) photoionization by the background field of distant active galactic nuclei and of 
local sources associated with the absorption systems. In this model, two gas phases with different 
ionization parameters T (ratio of hydrogen ionizing photon density to total hydrogen density njj) 
are commonly invoked in order to explain the presence of low and high ionization species, e.g., 
Mgn, Civ, and Ovi (Bergeron et al. |1994| ; Lu & Savage |1993j ). As T is defined by «h for a 
given ionizing field, one can also estimate the spatial extent of the ionized clouds along the line 
of sight (LOS) if the total hydrogen column density A r (H) is known. In general, the reliability 
of these models strongly depends on a good knowledge of the metal abundances involved. On 
the other hand, double LOSs toward background QSOs offer a unique possibility to resolve such 
clouds geometrically and determine transverse sizes. For instance, Smette et al. (1992) report 



lower limits of D = 0.7 to 2.2 hgQ 1 kpc for the diameters of two metal systems at z ~ 2 toward 
UM 673 AB, from the analysis of the line strengths in both QSO spectra. A somewhat different 
result, although with large uncertainties, is achieved through statistical simulations by Smette et 
al. (|1995| , hereafter Paper I), who infer 25 < D < 300 h^ 1 kpc for Civ clouds in the LOSs to HE 



1104-1805 AB. At higher resolution, Rauch ( 1997b ) has recently shown that density gradients on 
sub-kpc scales in gas associated with metal systems are not uncommon. These results suggest that 
C iv absorbers are composed of a large number of small cloudlets. 

HE 1104-1805 AB (z cm = 2.31, m B (A) = 16.7, m B (B) = 18.6, angular separation 8 = 3.195") 



was discovered in the course of the Hamburg/ESO Survey (Wisotzki et al. 1993), and has been 
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already studied at medium resolution by Smette et al. (Paper I). There is wide evidence for this 
QSO to be gravitationally lensed (Paper I; Wisotzki et al. |1995| ; Courbin, Lidman, & Magain |1998 ) 



so the proper separation between LOSs will depend on the (hitherto not well established) redshift 
of the lensing agent. Two papers aimed at detecting the lensing galaxy of HE 1104-1805 with 
somewhat different results have recently appeared. Using IR direct imaging observations, Courbin 
et al. ( 1998] ) estimate the lensing galaxy to be at z\ ens = 1.66, while Remy et al. ( |1998| ) find that 



their Hubble Space Telescope (HST) and ground-based direct imaging observations are consistent 
with z\ ens = 1.32. Neither of these results influences our photoionization models. Nevertheless, 
the small velocity differences and the similar line profiles between absorption lines in A and B at 
z = 1.66 (see § ||) probably exclude the damped Lya (DLA) system observed in A as the lensing 
agent, thus excluding the 2ii ens = 1.66 result. 

The primary aim of this study is to compare ionization conditions and metallicities in clouds 
at z a bs = 1-66 for the two LOSs toward HE1104-1805. We use new HST Faint Object Spectrograph 
(FOS), and ground-based New Technology Telescope (NTT) and Keck spectra. The Civ systems 
observed in the UV and optical spectra of A and B at redshifts close to zdla are very well suited 
to determine whether this gas is indeed photoionized: it can be assumed that they are associated 
with the high-density gas giving rise to the DLA system observed in A at z = 1.66162 and have 
therefore a common chemical history. Consequently, the relative element abundances should be the 
same in all these systems. Since the DLA gas is expected to be opaque to the ionizing radiation, 
it is possible to derive these element abundances without ionization corrections. 

To investigate the issue of ionization in the C iv systems, we have used the photoionization 
code CLOUDY (Ferland |i~993l) and the ionizing radiation field proposed by Haardt & Madau 



(1996). We have assumed that this radiation field also ionizes the gas that gives rise to the 
strong Ovi absorption observed along both LOSs. We have made a distinction between two 
gas phases with different ionization states: a low ionization phase, where both Mgn and Civ 
absorption occur, and a high ionization phase, where O vi but weak C IV absorption occurs. The 
photoionization models for each of these two gas phases are constrained by the observed column 
density ratios iV(Clv) to iV(Mgll) and iV(Ovi) to iV(Nv), respectively. 

Our paper is organized basically in two parts: estimation of the metal abundances in the 
DLA gas, and photoionization models for the Civ systems. The spectra are described in § ^. 
In § H| we describe the line fitting method used, and discuss the important role played by the 
continuum fitting to the HST spectra. § || is devoted to the metal abundances determined for 
the DLA system observed in QSO component A. § || presents possible physical scenarios for the 
C iv systems observed in A and B based on their association with the DLA system and on the 
line profiles. The CLOUDY models and the resulting physical parameters for the C iv and O vi 
absorbers are described in § ||. Finally, we outline our conclusions in § |7[ 
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2. Observations and Data Reduction 

An overview of the spectra used in this paper is displayed in Table [j]. We now detail the 
observations. 



2.1. HST Spectra 

UV spectra of HE 1104-1805 A and B were taken in November 1995 with the Faint Object 
Spectrograph onboard the Hubble Space Telescope. Target acquisition and spectroscopy were 
done using Grating G270H with the red detector and the 3."7 x 3."7 aperture. This configuration 
yields a spectral resolution of FWHM = 2 A and a wavelength coverage from 2222 A to 3277 A 
(Schneider et al. 1993j ). Total integration times of 1790 and 6690 seconds for QSO component A 



and B, respectively, resulted in variance weighted spectra of maximum signal-to-noise ratios S/N 
= 20 (A) and 17 per ~ 0.5 A pixel. 



2.2. Keck HIRES Spectra 

Optical spectra of HE 1104-1805 A and B were taken in January and February 1997 with the 
Keck High Resolution Spectrograph HIRES (Vogt et al. [1991 ) and an 0.86" slit at FWHM = 6.6 



km s . They range from 3620 to 6080 A. A full description of the observations and extraction 
method, as well as the spectra themselves will be presented elsewhere. In short, an image rotator 
was used to keep the slit off the second image, and as close as possible to the parallactic angle. 
The continua were matched by using the brighter A image continuum as a template, i.e., the 
continuum points of the B image spectrum were, in a sense, "stapled" to the A image continuum. 
This was done using polynomial fits in such a way that differences between the continua on scales 
larger than at average 300 km s _1 are divided out, but regions smaller than that retain their 
differences between the LOSs. The typical absorption line or absorption complexes between the 
spectra (which were omitted from the fit anyway) are not affected as can be seen from the strong 
differences in the metal lines despite the very similar Lya forest. We have also taken special care 
to omit the metal absorption line complexes from the continuum points, to be sure that we do not 
wipe out the differences. One-sigma arrays were derived from the Poissonian photon error, and 
rebinned onto the 0.04 A/pixel constant wavelength scale. 



2.3. NTT Spectra 

Optical spectra of HE 1104-1805 A and B were obtained in February 1996 with the echelle 
spectrograph of EMMI on the ESO New Technology Telescope in the Red Medium Dispersion 
mode under subarcsecond seeing conditions. The wavelength coverage is 3910 to 8290 A. Both 
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images were simultaneously centered in a 1" slit. Grating #9, grism 3 as cross-disperser, a F/5.2 
camera and a TEK 2048 2 CCD were used. This CCD provides 24^m pixels, corresponding to 0."27 
in the sky. The total integration time was 16 hours. The echelle orders were optimally extracted 
with a version of the extraction algorithm used in Paper I, modified to extract cross-dispersed 
spectra. The algorithm attempts to reduce the statistical noise in the extracted spectra to a 
minimum, and allows the correct separation of the two seeing profiles. It basically consists of 
the following steps: (i) a variance (Poisson statistics, read-out noise and cosmics) is assigned 
to each pixel. For each flat-fielded two-dimensional spectrum, two Gaussians of common width 
are simultaneously fitted to the profiles at each wavelength channel along the previously defined 
echelle orders using the Levenberg-Marquardt method (Press et al. 1986| ); (ii) the variation of the 
width and the position of the brighter component with respect to the orders in the dispersion 
direction is then fitted with low order polynomials; (iii) step (i) is repeated, this time with fixed 
width and position — given by the polynomial fits — thus allowing only the amplitudes to vary. The 
final, variance-weighted coadded spectra have FWHM = 0.8 to 1.0 A and maximum S/N ~ 77 (A) 
and ~ 27 (B). The continuum level was estimated separately in each order — skipping corrections 
for the blaze function — and for each QSO image. This was done by fitting low-order polynomials 
and cubic splines to featureless spectral regions. 



2.4. AAT Spectra 

Additional FWHM = 1.2 A resolution spectra of both QSO images taken with the 3.9m 
Anglo-Australian Telescope and covering the wavelength range 3170 to 7570 A were also used. 
They have already been presented in Paper I. 



2.5. Wavelength calibration of the HST spectra 

Special care has been taken to re-define the absolute zero point of the wavelength scale in 
the HST spectra. An off-center position of the targets in the aperture of the FOS may cause 
differences in the wavelength scale between the A and B spectra. By measuring wavelength 
positions of the galactic Mgn AA2796, 2803 lines, and assuming that this absorption takes place 
in the same cloud along both LOSs we found an offset of AX(A — B) = —0.42 A. This correction 
was applied to the B spectrum. Additionally, an overlapping 100 A wide spectral region around 
A = 3220 A in the HST and ^4^4T spectra allowed for a correction of the FOS wavelength scale 
to vacuum-heliocentric values. This was done by comparing the central wavelength differences of 
absorption lines in the A spectrum thought not to be blended. A relatively small correction of 
AX(AAT - HST) = +0.04 A was then applied to both HST spectra. 
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3. Absorption Line Analysis 

3.1. Continuum Definition in the HST Spectra 

Owing to the low FOS resolution, line blending introduces a serious problem when performing 
continuum fitting. This is particularly marked in the Lya forest, because of the lack of spectral 
regions free of absorption lines. For this reason, we decided to first determine the absorbed 
continuum — clearly dominated by two optically thin Lyman-limit systems (LLS) at z = 2.20 and 
2.30 and the DLA system (A) and a LLS at z = 1.66 (B) — before describing the intrinsic QSO 
continuum. Both spectra were corrected for galactic extinction (Seaton |197S| ) with E(B — V)=0.09 



(Reimers et al. 1995). The Lyman edges, arising from the superposition and blending of 
corresponding high-order Hi Lyman series lines, were modeled with Voigt profiles convolved 
with a FWHM = 2 A Gaussian describing the instrumental profile. The optical depth r at each 
Lyman break, given by the ratio of the extrapolated continuum to the Lyman one, determined 
iV(Hi), whereas b was mainly constrained by the shape of the edge. These parameters are listed 
in Table The total optical depth at each LLS is given by the Lyman continuum and the lines of 
the Lyman series. Note that the flux for A < 2200 A in the B spectrum is not completely absorbed 
(see Fig. p. 

The intrinsic QSO continuum of HE 1104-1805 A and B in the HST spectra for A obs < 3277 
A was found to be very well represented by two power laws (/ oc v a ) with a break at A = 2917 
A. To determine the power law parameters, a maximum likelihood fit was performed, matching 
simultaneously the intrinsic and the H I absorption continuum with the observed flux at selected 
absorption-free spectral regions including some regions dominated by significant Hi absorption 
lines. This led to the best-fit spectral indices a = —0.8, —1.5 for A, and —0.8, —0.9 for B (see 
Fig. |H). We believe this continuum estimation is a good one, since an extrapolation to longer 
wavelengths fits the scaled A AT spectra to well within la in regions of low Lya line density. 
Moreover, expected emission lines by Lya, Nv, and Ovi stand out well against the continuum 
as can be seen in Fig. [l|. Division of the flux by this continuum resulted in the normalized HST 
spectra shown in Fig. ||. 

The different continuum slopes of A and B, already pointed out in the discovery paper, might 
be a consequence of QSO component A being microlensed (Wisotzki et al. 1993| ). Furthermore, 



we find the spectrum of A bluewards of Lya emission to be softer than that of B, in agreement 
with the variability in the spectral slopes reported by Wisotzki et al. ( |1995| ) from the analysis of 
low resolution observations in the optical range made one and two years before ours. However, 
an alternative explanation for the different slopes might be differential reddening by dust grains 
in the DLA system observed in A, a possibility recently considered for the gravitationally lensed 



QSO 0957+561 (Zuo et al. |1997|) . 
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3.2. Line Profile Fitting 

In this section we describe the line fitting procedures used to obtain column densities 
of lines associated with the z = 1.66 absorption systems in A and B. In general, this is a 
nontrivial task because of the limited resolution of our spectra other than Keck HIRES, so that 
assumptions concerning line widths must be made. The lines were modelled with Voigt profiles. 
We distinguish between maximum likelihood fits, performed to lines in the Keck and NTT spectra, 
and "interactive" fits, performed to lines in the ^4^4Tand HST spectra (aside from Hi lines). 



3.2.1. Keck Spectra 

To determine column densities N and Doppler parameters b of lines in the Keck HIRES 
spectra, we x 2 -fitted Voigt profiles convolved with the instrumental profile to lines in the Keck 



spectra using the MIDAS program FITLYMAN (Fontana & Ballester |1995|) . These lines lie 
longward of Lya emission. Line parameters, i.e., rest-frame vacuum wavelengths, damping 
constants and oscillator strengths, were taken from Morton Q1991| ), and from Verner et al. (1994) 
for lines with revised /-values. 

Most of the column density errors (JiogN range from 0.01 to 0.10 dex; however, the smoothing 
introduced by rebinning the la-arrays sometimes underestimates the true flux uncertainties, 
making a\ og ]\r systematically too low. Thus, we made the following correction to the fit errors: 
for each line, we measured the amount of smoothing by calculating the ratio of the flux standard 
deviation from 41-pixel wide featureless regions in the data to the sigma from the error array. In 
the (few) cases where this ratio was larger than one, criogTV was corrected by this factor. Clearly, 
this applies only to unsaturated, resolved lines; a similar treatment to non-resolved complexes is 
less obvious because the <Ti og Ar's result from a more complicated Hessian matrix, and are no longer 
independent. We did not allow for this effect. 

Special attention must be given to the fits of lines associated with the DLA system observed 
in the spectrum of QSO component A because they will determine the metal abundances. To look 
for possible hidden saturation, we used the apparent optical depth method to compare apparent 
column densities -/V ap p(?;) of different transitions of a given ion in velocity space (e.g. eq. [1] in 
Lu et al. |1996j ; Savage & Sembach |1991| , for a description of the method). Although some lines 
are probably saturated (e.g. Alii A1670), in most of the cases where more than one transition is 
available we do not find significantly saturated structures, and the agreement between the fit and 
integration results is remarkable (cf. Table f|). 
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3.2.2. NTT Spectra 

In the NTT spectrum, the Mg n profiles at zdla are slightly asymmetric, suggesting that 
this system will also split into more components at higher resolution. Consequently, we fitted 
four-component Voigt profiles to these lines, with 6-parameters fixed at values found for the Sin 
lines in the Keck spectra. 



3.2.3. HST and A AT Spectra 

At even lower resolution, one cannot expect to recognize the line profiles properly; hence, 
to derive column densities of lines in the HST and A AT spectra, turbulence dominated line 
broadening was considered. Voigt profiles were interactively created and superimposed to the 
spectra using XVOIGT (Mar & Bailey |99|), while attempting to minimize the residuals. The 



redshift and Doppler widths used to create such line profiles were those found in a second fit with 
FITLYMAN of lines present both in the A AT and Keck spectra. Lines with asymmetric profiles 
were considered to be single, and the total column densities were fixed to the known Keck- values. 
In this fashion, the low and high-ionization species were distinguished using "low-resolution" 
Doppler parameters determined by the fits to Fen A1608 (b = 20 km s _1 ) and Civ A1548 (b = 44 
km s _1 ) lines in A, respectively. For B, the lines used were Alii A1670 (b = 20 km s _1 ) and Civ 
A1548 (b = 37 km s" 1 ). 

To estimate the uncertainties of our column densities we smoothed and rebinned lines 
observed in the NTT spectra to HST FOS resolution, and re-computed column densities with 
the procedure described above. The new column densities showed deviations of the order 0.1 
to 0.2 dex from the original, better determined values. Another source of error is our limited 
ability to de-blend metal lines from Lya forest lines. On the other hand, most column densities of 
transitions in the UV that contribute to the metal abundances are based on one line in the HST 
spectra and another in the A AT ones, e.g., Cn AA1036,1334; Oi AA988,1302 (see Fig. |). If these 
two effects compensate, we think that taking own = 0.2 dex for these ions is appropriate. 



3.2.4- Detection Limits 

We defined 3 a detection limits for metal lines in the HST and AAT spectra according to the 
formula (Caulet [1989]) 

FWHM , . 

where FWHM is the width of the spectral point spread function and <S/N> is the mean local 
signal to noise expected at the position of the line (measured as the inverse standard deviation 
of the normalized flux in small featureless stretches adjacent to the line). They range between 
W ohs = 0.13 and 1.00 A in the B spectra. 



- 9 - 



3.2.5. Hi Column Densities at z = 1.66 

To derive more accurate column densities for Hi, we used the normalized HST spectra. 
Because they completely cover the rest-frame spectral range down to 912 A for the z = 1.66 
systems, it is possible to measure iV(Hi) in both spectra more accurately than in previous 
studies on damped systems, by using higher Lyman series transitions. We simultaneously fitted 
two-component Voigt profiles to 11 resolved Hi lines in the normalized spectra of A and B using 
FITLYMAN. The fit solutions to lines in B were constrained by r = 2.34 at the Lyman edge. We 
estimate the neutral hydrogen column density of the DLA gas (spectrum A) to be log iV(Hl) = 
20.85±0.01; for the LLS (B) we obtained logiV(Hl) = 17.57±0.10. Notice that these values are 
independent of the ones estimated for placing the continuum and shown in Table II. 



4. The Damped Lya System Toward HE 1104 1805 A at z DLA = 1.66162 

Table |3| displays the fit results for lines associated with the DLA system toward HE 1104-1805 
A at zdla = 1.66162. A wide variety of singly and doubly ionized species is observed in this 
DLA system, but Civ and Siiv are also present. We now describe the Keck HIRES line profiles, 
referring to the left-hand panels of Figures |3|, |5] and |3], throughout this section. 



4.1. Low Ion Profiles 

The left-hand panel of Fig. |3| shows the line profiles of strongest low-ionization species in 
velocity space, relative to v = at z = 1.66164, the redshift of Mgn in the DLA system. These 
lines lie redward of the Lya forest. In each line complex, we have fitted four Voigt profiles with 
independent z, b and N values. The Nili and Fen results are based upon simultaneous fits to 3 
transitions; the Sin and Aim fits on 2 transitions; the Alii fit on 1 transition; and the Mgn fit on 
2 transitions (cf. Table |||). We fitted 4 components to the Mgn lines in the NTT spectra, with z 
and b tied to the values found for Si II. Only the three bluemost components show associated Zn n 



and Cm (see [4.3.1|) . 



From the high-resolution plots, we see that the low ions track each other quite closely, 
suggesting they occur in the same gas clouds. The whole profile is characterized by one cloud at 
v ~ —30 km s _1 with the strongest absorption, one cloud at v ~ +40 km s _1 with the smallest 
column densities, and two clouds with intermediate column density clouds lying in between. 
This "edge-leading asymmetry" seems to be a common feature of low-ionization absorption lines 
associated with damped Lya systems. It has been variously interpreted as a consequence of 
absorption by rotating gaseous disks (e.g., Wolfe et al. |1995b|) or as the signature of merging 



protogalactic clumps in hierarchical structure formation (Rauch et al. 1997a ). The different 



column density ratios at each velocity indicate clouds with different physical conditions (gas 
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density, metallicity, or even ionization) within 70 km s . 

4-. 1.1. Line Widths 

Remarkably, and despite the independent fits performed to each ion profile, we obtain fit 
solutions that uniquely characterize each cloud in redshift and broadening parameter (cf. Table ^ 
and Fig. ^). For instance, for the clouds at v ~ —30, —10, 10 and 40 km s~ x we find respective 
mean widths and standard deviations of (b) = 3.5 ± 0.7, 12.6 ± 2.7, 9.1 ± 2.2 and 7.6 ± 1.8 km s _1 
(fourth Nin component excluded due to large b- uncertainties; fit to Znii, Cm and Tin lines not 
included). Thus, we find a similar line-broadening mechanism in each of these absorption systems 
and, based on the small Doppler-parameter dispersions, cautiously favor turbulent gas motions as 
the dominant line-broadening mechanism. 

4.2. High Ion Profiles 

The left-hand panel of Fig. [| shows the Civ A1548, 1550, Siiv A1393, 1402 and, for comparison 
purposes, Alii A1670 velocity profiles associated with the DLA system. For the Civ absorption 
lines between [-80, +80] km s -1 , only four-component Voigt profile fits succeeded. Unfortunately, 
the poor S/N at the position of the Siiv lines and contamination by Lya forest lines do not allow 
a clear comparison with the Civ profiles. We decided not to fit these Siiv lines. Instead, we 
give upper limits for column densities based on the apparent optical depth method. The high ion 
profiles do not exactly track the low ion profiles: the data rather suggest at least part of the C iv 
absorption (velocity component 1) occurs in clouds without low ionization species. On the other 
hand, C iv velocity components 2 to 4 show a certain resemblance to the edge-leading asymmetry 
of the low ion profiles; however, our fit solution yields relatively large 6-values, suggesting 
that C iv, regardless of which line broadening mechanism dominates — thermal or turbulent gas 
motion — occurs in hotter gas regions than the low ion clouds. 

4.3. Abundances 

Metal abundances in the DLA system, normalized to solar values and defined by 

[M/H] = log(JV(M)/JV(H)) - log(JV(M)/JV(H)) (2) 

were computed using the column densities integrated between [-50,60] km s _1 . They are listed in 
Table |j. Solar abundances were taken from Verner et al. ( |1994| ). Singly ionized species provide the 
bulk of the total element column densities, except for O I, assumed to be the dominant ionization 
stage of oxygen. 
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4-3.1. Zn, Cr and Ti Abundances and Dust 

Fig. U shows the Keck velocity profiles of the most outstanding Znii and Cm transitions. 
The fit solutions lead to three-component profiles for the Cm complex at v = —29.1, —7.3 and 
7.4 km s , and two-component profiles for Znii at v = —30.0 and —9.0 km s _1 , relative to 
z = 1.66164. The column density ratios relative to solar vary from iV(Znii)/iV(Cr Ii) = 3.47 
(bluemost component) to 1.73. If the Zn and Cr abundances as well as the ionization level were 
the same in these clouds (which is very probable), then this variation would indicate that the 
dust-to-gas ratio in this DLA gas is inhomogeneous within 20 km s _1 , with a higher dust content 
in the higher density Znii component. 

The variation of the abundance ratios of refractory elements among different clouds associated 
with the DLA gas provide insights into the presence of dust in the disk of damped Lya galaxies 
(Lu et al. |1996|), because dust grains can be locally destroyed by passage of supernova shocks 



(Sembach & Savage 1996, and references therein). We can extend our analysis of elemental 
abundance ratios to iron and nickel, also expected to be depleted into dust. The column density 
ratios Znii to Fell and Znii to Nin relative to solar are respectively Zn/Fe = 5.5 and Zn/Ni 
= 7.9 for velocity component 1, and Zn/Fe = 2.5 and Zn/Ni = 2.8 for velocity component 2, 
thus in concordance with what one observes for Cr in the corresponding clouds (we estimate the 
corresponding uncertainties to be no larger than 0.4). Although these variations are small, 0.3 - 
0.4 dex, it seems that in the v ~ —30 km s _1 cloud the effect of dust depletion is more important 
than in the cloud at v ~ —10 km s _1 . On the other hand, as pointed out in section [4.1.1| , 
the bluemost component is characterized by narrower lines than component 2 in all the ions 
considered (including Znii and Cm). It would be of great interest to discern whether such line 
width differences have a thermal origin (contrary to what was stated in |4.1.1| , however), because it 
would give evidence for dust depletion being more effective in cooler gas. 



Lu et al. ( 1996 ) have presented arguments for a pure nucleosynthetic origin of the elemental 
abundance pattern observed in Lya damped systems at low metallicity. These arguments are: 
the low N/O ratio, which we also find in this DLA system (see next section); the a-element 
overabundance relative to Fe-peak elements, also observed for the abundance ratios of Si to Fe, 
Cr, Mn and Ni in our Keck HIRES data; and the underabundance of Al relative to Si and Mn 



relative to Fe (the "odd-even effect"; see Lu et al. 1996 for details), which we do not find in this 
DLA system. Instead, we derive [Al/Si]= +0.25 ± 0.15, and [Mn/Fe]= -0.02 ± 0.14 (cf. Table g), 
although we recall that the Al abundance is only based on the Al n A1670 line. Given the high 
Mn/Fe ratio and the argument given in the last paragraph, we suggest that stellar nucleosynthesis 
alone is not likely to produce the relative abundance pattern observed in this DLA gas. 

Based on the total zinc abundance [Zn/H] = —1.02 ± 0.01 relative to the solar value, we 
derive a metallicity Zdla — 1/10 Zq for this system. This zinc abundance is somewhat lower than 



the value [Zn/H] = —0.8 reported by Pettini et al. ( 1997 ) [same data as in Paper I], who also 



used log(Zn/H)Q = —7.35, and more accurate since it is based on a simultaneous fit to two Znii 
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lines.^ Additionally, we deduce from [Cr/H] = —1.46 ±0.02 a dust-to-gas ratio of ~ 0.11, using 
the definition given by Vladilo ( |1998[ ), eq. [19], and considering that most of chromium should 
be incorporated into dust grains in the ISM of galaxies showing damped Hi absorption (e.g. 
Pettini et al. 1990| ). Furthermore, titanium, another refractory element, is found to have [Ti/H] 
= —1.50 ± 0.07, based upon the unsaturated Tin AA1910.6, 1910.9 lines. This is in full agreement 
with the incorporation of Ti and Cr into a dust phase in the neutral ISM of this damped Lya 
galaxy. Thus, based on these abundances, we conclude that (1) this DLA system does not differ 
too much from other ones at higher redshifts, given the scatter observed in [Zn/H] (cf. Fig. 3 in 
Pettini et al. 1997|) ; (2) there is evidence for the presence of dust. 



4.3.2. The Abundance Ratio O/N 
The crucial abundance ratios in this work are [O/N] and [C/Mg]. 

Based on the Oi AA988,1302 lines, we find [O/H] = —0.98 ±0.20, in very good agreement with 
the abundance ratio of Zn. This supports O I as a very good tracer of H 1 — as can be expected 
from its ionization potential, 13.62 eV — and suggests that O is not depleted into dust in the ISM 



of this damped Lya galaxy, in agreement with observations in the local ISM (Cardelli et al. 1991) 
Moreover, since both O and Si have been observed to have solar abundance ratios in Galactic 
halo stars and in metal-poor dwarf galaxies, one would expect [Si/O] ~ in DLA systems (Lu et 



al. 1996 and references therein). We do obtain the same abundance ratios for O and Si within the 
errors, and [Si/H] is based on reliable column-density measurements of two Sill lines (observed 
in the Keck HIRES spectrum), making our oxygen abundance estimation yet more confident. 
Concerning the abundance ratio of nitrogen, the bulk of [N/H] is provided by the column density 
of Nil A1083. This line is very probably not contaminated by a Lya forest interloper, given the 
absence of absorption at the same wavelength in the spectrum of B (see Fig. ||). In addition, 
nitrogen — like oxygen — is also expected not to be depleted in the ISM (Cardelli et al. 1991| ). In 



consequence, we are confident of an abundance ratio of [0/N]= 0.88 in this DLA gas, that is, 
O/N is 8 times greater than the solar ratio. This is in qualitative agreement with observations of 
damped Lya galaxies at higher redshifts (Pettini et al. 1995; Lu et al. 1996j ), and with galactic 



chemical evolution models, because these two elements have different nucleosynthetic origins, 
oxygen being produced in much shorter timescales than nitrogen. 



4 However, our 40 km s 1 resolution NTT spectra yield Zn and Cr abundances completely consistent with the 
Keck results, thus validating abundance studies of Zn at medium resolution. 
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4-3.3. The Abundance Ratio C/Mg 

The carbon abundance is based on the fits to the Cn AA1036, 1334 lines. C/Mg is found to 
have the solar value within the errors, but iV(Mg) might be underestimated through saturation of 
the Mgn lines in A. 



5. Geometry of the Absorbers at z = 1.66 

Fig. | shows the velocity profiles of the Civ A1548,1550 and Mgn A2796,2803 doublets (at 7 
and 40 km s^ 1 resolution, respectively) toward HE1104-1805 A (left) and B. In both panels, v = 
km s _1 corresponds to z = 1.66164, which is the redshift of Mgll in the DLA system. We have 
arbitrarily numbered the Civ complexes at z = 1.66143 (1), z = 1.66280 (2) and z = 1.66465 (3) 
in the A spectra, and at z = 1.66184 (4), z = 1.66284 (5) and z = 1.66493 (6) in the B spectra, so 
they will be refered to as systems 1 to 6 throughout the following sections. Also shown in Fig. ||, 
covering a larger range of velocities, are the profiles of the Lya, Ly/3 and Ly7 Hi lines. Note that 
Hi presents (at least) two components both in A and B. 

C IV systems 1 to 4 show associated Mg II, although shifted in velocity, as is more evident in 
systems 1 and 4 (the line shapes suggest that either of these systems will probably split into more 
components at even higher resolution). However, most of the Mgn seen in A at v = is due to 
the DLA system. In B, system 4 is identified with the LLS. The presence of Mgll in absorption 
systems 2 and 6 is less evident at this S/N, so only upper limits can be derived. Civ system 5 will 
not be considered here. 

Do C iv and Mg n occur in the same clouds? Although the present data show a correspondence 
in velocity, there might not be a physical association between both ions. However, we will show in 
the next section that photoionization simulations of these clouds do indeed predict the presence 
of both ions in a common gas phase, in agreement with our observations. Furthermore, we know 
from studies at high resolution that line profiles of low and high ions do track one another in 
Lyman- limit systems (Prochaska & Wolfe 1996j ). As discussed in § 6.2 it is also possible that part 



of the C iv arises in the same highly ionized gas that gives rise to O vi. In any case, the velocity 
profiles of system 4 (LLS) show that, if a physical association of both ions is correct, at least part 
of the Civ arises in clouds where no Mgll is present (see also the right-hand panel of Fig. ||). 

From inspection of Fig. ||, it seems likely that LOSs A and B cross common absorbers. This 
is suggested by the similar line profile pattern in both spectra. If a common absorption complex 
gives rise to systems 1 (A) and 4 (B), then the different equivalent widths of the Mgn lines in 
A and B, in contrast to the more similar C iv equivalent widths, suggest gas inhomogeneities on 
spatial scales smaller than the linear separation between LOSs S± = 8.3 h^Q 1 kpc, for qo = 0.5 and 
a lens at z = 1, thus suggesting that Civ arises in a more extended region than Mgll. However, 
we find that the velocity difference between the Civ clouds in A and B is (Av) = —11 ± 2 km s _1 
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while (Aw) = — 1 ± 1 km s _1 for the corresponding Mgn clouds. If these velocity differences are 
a consequence of peculiar cloud motions, then a physical association of C iv and Mg n is still 
compatible with the data. Since it is not the aim of this study to determine the origin of both this 
velocity difference and the ~360 km s _1 velocity span of the Civ absorbers along both LOSs, we 
have analyzed each system separately. 

There are two alternative interpretations for the Hi column densities in A and B: (1) the DLA 
system arises in a disk-type galaxy (Wolfe |l995a|) with the LOS to A passing through the gas in 
the halo and the disk and LOS B passing through the halo gas only, or (2) in models of hierarchical 
structure formation the DLA system arises in the central region of a protogalactic clump and 
LOS B crosses the surrounding less dense gas at an impact parameter of a few kpc (Rauch et 



al. 1997a ). Our data do not allow us to discriminate between these two models. Consequently, 
in the following we will simply consider C IV absorption systems 2 to 6 to arise in clouds in the 
extended halo of the cloud giving rise to the DLA system (system 1), giving explicit references to 
one of these models. 



6. Ionization State and Chemical Composition 



We have used the photoionization code CLOUDY (version 84.12a; Ferland 1993| ) to investigate 



the ionization state and metallicity in the LLS observed in the B spectra of HE 1104-1805 
(system number 4 in Fig. |^) and in three further Civ-Mgn systems observed in A (systems 2, 
3) and B (system 6). The clouds giving rise to these systems are represented by parallel slabs 
illuminated on one side by the radiation field J u proposed by Haardt & Madau ( |1996| ), consisting 
of the background flux contributed by QSOs and AGNs, which is attenuated by H and He 
absorption in Lyman- limit systems and Lya forest clouds. This radiation field also has a diffuse 
component due to recombination continuum radiation. Our model assumes J(912) = 0.37 x 10~ 21 
erg s _1 cnr 2 Hz -1 sr _1 at z = 1.66. Although this geometry does not perfectly describe a cloud 
being illuminated by an isotropic incident radiation field, it should not introduce errors larger than 
a factor of ~ 2 (Bergeron and Stasinska 1986). In particular, cloud sizes along the LOSs resulting 



from this model need to be considered upper limits, because considering slabs illuminated only on 
one side underestimates the true ionizing radiation field. 

As we argue below, the wide variety of low- and high-ionization stages present in these 
systems makes it necessary to model the gas clouds with two zones of different ionization levels: 
(1) a "low-ionization phase", where absorption by singly, doubly, but also triply ionized atoms 
occurs^; (2) a "high-ionization phase", where O vi absorption occurs. However, we will also discuss 
the possible existence of a third, "intermediate-ionization phase" . 

The CLOUDY simulations predict column densities of the expected ionization stages, suitable 



5 The term "low" has specifically the purpose to distinguish both gas phases. 
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to be compared with the observations. As a general strategy, we attempt to reproduce the 
observed ionic column density ratios iV(Crv) to iV(Mgii) in systems 2, 3, 4 and 6, and the ratio 
iV(Ovi) to iV(Nv) in the high-ionization phase by varying the ionization parameter T. We assume 
that both gas zones are ionized by the Haardt & Madau radiation field, and have the same relative 
abundances as found in the DLA gas. Since the column density ratios are quite insensitive to the 
metallicity of the gas [M/H] for a wide range of ionization parameters, we can constrain [M/H] 
using the individual observed column densities. 

Clearly, the assumption of same relative abundances in A and B must not hold for gas-phase 
abundances of elements known to be depleted by condensation into dust grains. If the dust-to- gas 
ratio in the halo of this damped Lya galaxy were considerably lower than in its DLA region (or 
even zero), photoionization models should arrive at ionic column densities of refractory elements 
that are underestimated, compared with the observed value. We observe this effect for Fen (see 
next section). 

6.1. Low-ionization Phase 

6.1.1. The Lyman-limit System at z = 1.66184 toward HE 1104-1805 B 

Table |5| displays the observed column densities of ions associated with the LLS at z = 1.66184. 
The low and high ion velocity profiles are shown in Figures || and respectively (right-hand 
panels). Besides the DLA system, this is the second most metal-rich system. Prominent ions 
observed are: Mgn, Alii, Aim, Sin, Sim, Cn, Nm and Fen. The detection of Fen is marginal 
but significant. Due to the low S/N and extraction artifacts at the position of Fen A1608 (Fig. ||, 
top right spectrum), we decided to fit one-component Voigt profiles to the three strongest Fen 
lines in the NTT spectra. The detection of Al n and Al III is qualitatively consistent with the 
assumption of DLA relative metal abundances in the LLS. 

Taking the arguments given in § || into account, in order to perform photoionization 
simulations for the LLS we have to first determine the range of possible values within which 
A?"(Crv)/iV(MgIl) is likely to vary. Since the line profiles of Mgn do not trace the whole velocity 
range of Civ, a very conservative upper limit of ~ 20 for the ratios is given by the total column 
densities of the fitted Voigt profiles. However, from the right-hand panel of Fig. ||| we see that the 
low ion profiles (here represented by Al n A1670) coincide quite well in velocity space with C IV fit 
components 2 and 3. Considering that these Mgn line profiles at Keck HIRES resolution would 
not differ too much from the Alii A1670 ones, and taking the column densities integrated over 
these components leads to iV(Crv)/iV(Mgii)= 5.3, a much more realistic column density ratio. 

Table || displays the column densities predicted by CLOUDY for the LLS using the Haardt 
& Madau ionizing radiation field (column labelled MODEL 1). Assuming for this system DLA 
gas-phase abundances, a column density ratio iV(Clv)/iV(Mgll)= 5.3 implies lognn = —1.8 cm -3 



-16- 



(or logT = —2.95 for this radiation field); in other words, the gas is relatively highly ionized 
with iV(Hll)/JV(Hl)w 200, but Mgn is still present. The observed logiV(Hl) = 17.57 leads to a 
typical (model dependent) cloud size along the LOSs of Su < 1.6 kpc. For comparison purposes, 
a power law of the form / oc v~ 2 was also used as ionizing background (column labeled MODEL 
2 in Table |6|). A harder radiation field is not able to reproduce A r (Civ)/A r (Mgii)= 5.3 with the 
assumed element abundances. However, even the selected power law model fails to simultaneously 
reproduce C n and C iv, while the Haardt & Madau model does, due to the continuum break at 
the Hen edge. 

Our model assumes that the LLS gas is in photoionization equilibrium. Deviations from this 
equilibrium can lead to dramatic underestimations of the cloud lengths parallel to the LOSs by 
photoionization simulations, because the neutral hydrogen fraction is overestimated (Haehnelt et 
al. |1996a ). However, at the density derived for this LLS, bh ~ 0.01 cm -3 , significant departures 



from the equilibrium temperature (T = 1.5 x 10 4 K) are rare, as is shown by Smoothed Particle 
Hydrodynamics simulations (Rauch et al. 1997a| , Haehnelt et al. 1996b). Thus, the hydrogen 



recombination timescale in this regime, ~ 6 Myr, is short enough to allow line cooling to balance 
photoheating processes. Consequently, we believe the CLOUDY sizes derived for this LLS to be 
reliable. 



Gas Metallicity in the Lyman-limit System. The photoionization models described above 
are quite independent of the gas metallicity [M/H] over the whole range of possible gas densities; 
therefore, [M/H] can be determined by matching predicted and observed column densities. We 
find that, regardless of which model is assumed, [M/H]lls = [M/H] DLA — 0.2 represents the best 
prediction for nine ions observed in this system (Alii, Aim and Fen are not considered). In 
particular, A^(Cn), A(Civ), A(Mgn), and A(Sim) are simultaneously very well reproduced if 
Zlls = 0.63 Zdla- This is an upper limit for Zlls because due to saturation of the Mgn lines 
in A, we obtain only a lower limit for [Mg/H]; hence, to reproduce A r (Mgii) in B, a larger Mg 
content would require an even lower metallicity in the LLS relative to the DLA gas. Varying the 
relative abundances within the observational errors in the column densities leads us to estimate 
that this result is significant at the la level (for comparison, considering solar relative abundances 
in the LLS leads to [M/H] LLS = -1.5). 

Studies of our galaxy halo gas show no systematic differences in the gas-phase abundances 
within galactocentric distances of 7 to 10 kpc in various directions, suggesting uniform physical 
properties over these radii. In addition, warm disk clouds show gas abundances that are 0.2 to 0.6 



dex lower than those in warm halo clouds (Sembach & Savage 1996 , and references therein; Savage 



& Sembach 1996 , but see Cardelli et al. |1995[) . If we assume that the LLS observed at z = 1.66184 



in HE 1104-1805 B arises in the halo of the z = 1.66162 damped Lya galaxy seen in A, a negative 
gradient in metallicity from DLA to halo gas implies that this halo gas has not yet been fully 
enriched with metals. This might be a consequence of different star-formation rates, in which case 
LOS B would be probing gas regions of lower star-formation rate than LOS A. Alternatively, such 
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an abundance gradient could also be explained by a disk-to-halo gas (and dust) transfer yet in 
early stages, if gas in the halo originates in the disk. 



Fe and Al Abundances in the Lyman-limit System. Singly ionized iron is the most 
outstanding outlier in our model, yielding an Fen abundance one order of magnitude lower 
than observed (cf. Table |6|). We conclude that this difference can only be due to different 
iron gas abundances between A and B, qualitatively consistent with the absence of dust in 
the halo of this damped Lya galaxy. This situation resembles that in our Galaxy, where the 
degree of dust depletion in halo clouds is smaller than in Disk clouds (Sembach & Savage 1996 ). 
Further CLOUDY simulations show that the observed log iV(Fe n)= 12.47 can be reproduced if 
[Fe/H] LL s = -0.8 (Haardt & Madau |1996|) or -1.0 (f u cx v~ 2 ). An opposite effect is found for 
aluminum, where our model reproduces Aim only if [Al/H]= —1.0, that is, if the Al overabundance 
is limited to the DLA gas (but such overabundance might be a consequence of a saturated Alii 
A 1670 line in A [see also section |3.2.1| ). 



A third ionization phase? N(C iv)/iV(Mg n)= 5.3 also requires A^Cin) to be lower by 1.5 a 
than observed. This suggests that there might in fact be a third "intermediate-ionization gas 
phase" , where part of the C iv, C in, N in and Si iv but no singly ionized species occur (neither 
does Si in, whose ionizing potential of 33.5 eV is considerably lower than that of Cm). We would 
be thus observing blends of lines arising in two phases at similar redshifts. The existence of such 
a gas phase is fully consistent with the model predictions for Si III and Al III in the low-ionization 
phase (provided the Al relative abundance is lower than in the DLA gas) , and with the C iv line 
profiles, showing a much wider velocity span than the low ions (Figg). Unfortunately, the HST 
spectral resolution does not allow an appropriate analysis of the Cm and Nm line profiles. 



6.1.2. The Civ Systems at z = 1.66280 and z = 1.66465 toward HE 1104-1805 A 

The iV(Civ)/iV(Mgii) ratios in the Civ systems at z = 1.66280 and z = 1.66460 in A 
(systems 2 and 3 in Fig. ^) are relatively well constrained by the observations, so we basically 
repeated the procedure described in the previous section, i.e., we searched for CLOUDY solutions 
that reproduce these ratios by varying T and Z. Since only the total neutral hydrogen column 
density is known, iV(Hi) was distributed according to the 7V(Mgn) ratios. 

Table ^ displays the column densities predicted by CLOUDY for System 3. Again assuming 
DLA gas-phase abundances, and the Haardt & Madau metagalactic radiation field, we find for 
this system that the observed A^(C lv)/iV(Mg n)= 70.8 can be well reproduced if lognn = —2.21 
cm -3 (or logT = —2.54 for this radiation field), implying a longitudinal cloud size of S\\ < 0.9 kpc, 
where the inequality stands for [Mg/H] > —0.97. This size estimate assumes both Civ and Mgn 
to arise from the same density region in photoionization equilibrium. 



-18- 



In System 2, the detection of Mgn is uncertain, and we can only derive 
iV(Clv)/iV(MgIl)> 117.5, which requires that logn H < -2.30 (logT > -2.45), or S\\ £ 
0.5 kpc for iV(Hl)= 16.05. Predicted column densities for system 2 are shown in Table f^. The 
disagreement for iV(C in) arises as a consequence of a bad fit to the C in A977 line. 

None of these photoionization models necessarily requires that Z < Z-qla- 



6.1.3. The Civ Systems at z = 1.66493 toward HE 1104-1805 B 

Observed and predicted column densities for the Civ system at z = 1.66493 in B (system 6) 
are shown in Table |9|. Because the detection of Mg n is very uncertain, our model was required 
to reproduce Civ assuming Z = 0.63 ^dla- We arrived at lognn = —2.40 (or logT = —2.35), 
implying S\\ ~ 0.6 kpc. 



6.2. Ovi-phase 

Possibly the most interesting result of this study is the significant detection of strong O VI 
absorption at z = 1.66 in both LOSs toward HE 1104-1805. Fig. |?] shows the corresponding section 
of the HST spectra of the A and B images with tickmarks above the spectrum indicating the O VI 
AA1031,1037 doublet at z = 1.66253 and the Cn A1036 lines associated with the damped Lya and 
Lyman-limit systems. Tickmarks below the flux level indicate lines identified with other systems. 
The O VI A1031 lines at A = 2747 A have similar rest-frame equivalent widths of W T = 0.93 db 0.07 
A (A) and 0.86 ± 0.10 A. Besides Sim A1206 at z = 1.28, neither further metal lines nor a Ly/3 
line has been identified at this wavelength, but contamination by a weak Lya line is not ruled 
out.0 Voigt profiles convolved with the instrumental profile have been overplotted in the spectra 
of A and B. They have Doppler parameters b = 180 km s -1 (A) and 110 km s _1 (B) (obtained 
from Gaussian profile fits to the Ovi A1031 lines) and common column density log A?" = 14.95. N 
shows little variation with b within 110-180 km s _1 in this region of the curve of growth. Clearly, 
these large Doppler values do not necessarily represent the true line widths. Rather, the observed 
line profiles are probably made up of more than one O vi line (but see next paragraphs) . For 
comparison, the dotted lines show the resulting profiles of 3 (A) and 2 (B) O vi doublets with 
common b = 50 km s _1 and total column density log N = 15.1 (A) and 15.0 (B). The doublets are 
placed at redshifts corresponding to systems 1, 2, and 3 in A, and 4 and 5 in B (cf. Fig. ||). 

Although Ovi has been shown to be relatively common in QSO absorption systems (Bahcall 



et al. 1993; Bergeron et al. 1994; Buries & Tytler 1996j), its detection along the two LOSs to 



6 It seems very unlikely that the absorption feature at A ~ 2747 A is entirely due to a Lya forest interloper. In 
such a case, two Lya lines with N(Hl) ~ 10 16 cm~ 2 (b = 30 km s _1 ) would be required to reproduce the absorption 
profiles; however, no associated C IV is found, as it would be expected in this class of Lya clouds. 
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HE 1104-1805 at z = 1.66 allows us to directly prove for the first time that this ion does indeed 
arise in very extended gas clouds. There were indications of this for systems at lower redshift 
(e.g. Bergeron et al. 1994 ; Lu Sz Savage 1993| ) but no firm evidence. Large, low-density Ovi 



gas clouds confirms the prediction by Rauch et al. ( 1997a| ) based on simulations of protogalactic 
clumps (PGCs). In this model, PGCs with a few times 10 9 M0 (in baryons) are embedded in a 
filamentary, low density and largely featureless Ovi phase with spatial extent of up to several 
hundreds kpc at the iV(Ovi) = 10 13 cm" 2 column density contour (cf. Fig. 3 in Rauch et al.). 
Detectable O vi is indeed expected to be more extended than C iv. The gas is photoionized, and 
the lines are bulk motion broadened and wider than C iv, as they sample the peculiar velocities 
over a larger volume. 

Observational evidence for the physical processes giving rise to Ovi — collisional or 
photoionization — remain, however, debatable so far, basically because contamination by Lya lines 
makes it difficult to resolve the Ovi AA1032,1036 doublet profiles appropriately. In our case, a 
quantitative assessment of the nature of this gas phase is possible, since no N v is detected at this 
redshift (Fig. |8|), and photoionization models can be well constrained under certain assumptions. 

At the hydrogen density nn ~ 0.01 cm -3 derived for the ClV+Mgil clouds, the fraction of 
ionizing photons with the E > 114 eV necessary to ionize enough Ov into "observable" Ovi is 
too small. Consequently, an additional, less dense phase is required to explain the strong Ovi 
absorption observed in A and B. Since a unique interpretation of the line profiles is difficult 
with the present HST data, we made the following, simplifying assumptions: (1) both LOSs 
pass through a common phase giving rise to the Ovi absorption lines, suggested by the similar 
equivalent widths and the zero velocity difference between lines in A and B; and (2) the absorption 
lines arise in one cloud, suggested by the line profile symmetry (this assumption may not be valid, 
but it is equivalent to integrating iV(Ovi) over several line components). A 3a upper limit of 
W T = 0.034 A for the Nv A1238.8210 line (the stronger line of the N v doublet) measured in the 
AAT spectrum of A leads to a 3a upper limit of log A^(N v) < 13.2. The non-detection of N V 
places the stringent (but otherwise very conservative) lower limit of N(0 vi)/iV(N v) > 60. This 
ratio requires densities lower than log ran = —3.55, or a large ionization parameter logT > —1.2, 
if we assume the same relative abundances as found in the DLA gas. This is a striking result, 
given that a lower (closer to solar) [O/N] ratio would lead to even less dense Ovi clouds. In this 
regard, we must bear in mind that the CLOUDY simulations above depend partly on the relative 
abundances assumed. For instance, variations of the [O/N] ratio within the estimated column 
density uncertainties will lead to differences in the predicted parameters. 

It is nontrivial to obtain an overall picture of this phase because the total hydrogen column 
density is not known. As a consequence, photoionization models will not reproduce the individual 
observed column densities iV(Ovi) and JV(Nv) uniquely, but both iV(H) and the gas metallicity 
Z will determine them. This is shown in Fig. [9|, where we have plotted all the pairs (Z,N(H)) 
that reproduce log iV(0 vi) = 15.0 and log JV(Nv) = 13.2. We can see that, in general, a lower 
metallicity of the highly ionized gas relative to that of the DLA gas, say Z < 0.63Zdla 5 will 



-20- 



be compatible with the observations if iV(H) is large, logiV(H) > 20.0. Additionally, we find 
that, regardless of the gas metallicity, the clouds where this phase occurs must be at least one 
to two orders of magnitude larger than the systems giving rise to Mg n absorption. These cloud 
sizes (5|| >~ 100 kpc), if representative of transverse dimensions, are widely consistent with the 
detection of O vi of similar strength in both LOSs. Typical masses derived from these models are 
~ 10 9 M . 

To see if Civ is expected in such a model, we performed further ionization simulations for 
larger ionization parameters, requiring that the assumed total hydrogen column density and 
metallicity always led to logiV(Ovi) = 15.0, i.e., curves departing from points on the Z vs. iV(H) 
curve in Fig. ^. The predicted N(C iv) column densities in such models are shown in Fig. |l0| for 
the Haardt & Madau (solid line) and a power-law radiation field (dashed line). At the minimum 
value of r allowed by the observations logT = —1.2, and for a wide range of gas metallicities, 
CLOUDY predicts log./V(Civ)« 14.4, in agreement with the observed value. However, at higher 
ionization parameters, e.g., A^(0 vi)/iV(N v) ~ 90 or logT = —1.0 (corresponding to a 2a 
significant non-detection of Nv), iV(Clv) becomes ~ 40% smaller, independently of Z. This 
result in turn indicates that at least part of the observed C IV arises in highly ionized gas, with the 
remaining contribution coming from the low-ionization gas phase. The present HST data do not 
enable us to quantitatively assess this issue and we can only conclude that C iv is likely to occur 
in both phases. 

Collisional ionization of gas in thermal equilibrium can also reproduce the observed 
N(Ovi)/N(Nv) ratio at T = 10 5 ' 45 K if one assumes solar relative abundances (Sutherland & 
Dopita p93| ), so we cannot rule out this process as responsible for ionizing Ov into Ovi. In 
such a case, the line broadening would be mostly due to macro-turbulence. Collisionally ionized 
O VI gas in a Lyman-limit system at z ~ 3.4 has been recently proposed by Kirkman & Tytler 



(1997), who detect both Civ and O vi at the same redshift in high resolution data. These authors 
favor collisional ionization due to the high kinetic temperatures implied by the line widths of both 
ions. In the case of HE 1104-1805 AB, additional, higher resolution (R ~ 23 000) UV spectra are 
needed, in order to finally discriminate between collisional and photoionization as the dominant 
ionization mechanism in the Ovi-phase. 



7. Conclusions and Final Remarks 



We outline our conclusions as follows: 



1. The damped Lya system observed at z = 1.66162 in the ultraviolet and optical spectra 
of HE 1104-1805 A has a neutral hydrogen column density of logiV(Hi)= 20.85 ± 0.01. 
The 6.6 km s" 1 resolution line profiles show that this high-density gas is distributed 
in four clouds spanning ~ 80 km s . Since the low ionization and Aim line profiles 
are observed to track fairly well, these ions are likely to arise in common clouds. Their 
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b- values are consistent if the broadening mechanism is not purely thermal. The C iv 
line profiles show that at least some of Civ absorption occurs in clouds with no singly 
ionized species. We find a metallicity of Zdla — 1/10 Z & and a dust-to-gas ratio 
of ~ 0.11, based on [Zn/H]= —1.02 and [Cr/H]= —1.46. Other element abundances 
are: [C/H]= -0.91, [N/H]= -1.86, [0/H]= -0.98, [Mg/H]= -0.97, [Al/H]= -0.77, 
[Si/H]= -1.02, [Ti/H]= -1.50, [Mn/H]= -1.61, [Fe/H]= -1.59 and [Ni/H]= -1.68. 
The O/N ratio, 8 times larger than the solar value, is consistent with galactic chemical 
evolution models. However, the observed elemental abundance pattern must be modified by 
condensation into dust grains given (a) the underabundance of Cr and Ti relative to Zn; 
and (b) the variation of the abundance ratios of these and other refractory elements among 
different clouds associated with the DLA system. 



2. Photoionization by the metagalactic radiation field proposed by Haardt & Madau ( |1996 ) 



implies that the Lyman-limit System observed at z = 1.66184 in the spectra of HE 1104-1805 
B and the three further C iv systems observed at similar redshifts in A and B all belong to 
the same category of absorbers, namely small (LOS sizes 0.5 to 1.6 kpc) and relatively dense 
(tih = 0.01 cm" 3 ) clouds where both Civ and Mgll absorption occurs. 

3. We find evidence for the LLS and one further Civ absorption system observed in B to have 
lower metallicities than observed in the DLA gas, Z = 0.63 Zdla- This result is significant 
at the 2 a confidence level if the column density uncertainties are not underestimated. 

4. In the context of galaxy formation, the observed Hi column densities suggest that LOS B 
intercepts gas in the "halo" of a protogalaxy at z = 1.66, while LOS A crosses its denser, 
central gas regions. If this picture is correct, then we are observing metal-poor gas in the 
halo of such a galaxy, compared to the high-density gas (within a transverse separation of 
~8 h^Q kpc). This in turn suggests regions of different star-formation rates or, alternatively, 
metal enrichment through gas transfer from the inner to the outer regions of the protogalaxy 
still in early stages. 

5. The presence of a highly ionized gas phase giving rise to the observed Ovi absorption of 
similar strength in both spectra is evident. If photoionization is assumed, the observed 
log iV(0 vi) = 15.0, and the non-detection of Nv implies large (~ 100-200 kpc), low-density 
(~10 cm -3 ) clouds where strong Ovi but weak Civ absorption occurs. On the basis of 
these results, we suggest that this highly ionized gas surrounds the clouds giving rise to the 
damped Lya and Lyman-limit systems. Such an extended low-density Ovi phase confirms 
the predictions of simulations of protogalactic clumps (Rauch et al. 1997a). 



Concerning the negative gradient in metallicity from gas crossed by LOS A to gas crossed 
by LOS B, it must be pointed out that such an effect, if real, does not necessarily support the 
disk/halo scenario. Indeed, it is still not clear whether present-day rotating disk- like galaxies 
should have had a similar appearance in their early stages of formation, nor is it well understood 
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what we should define as protogalactic "halos" . A good deal of progress is being made to resolve 
this paradigm: a link between metallicity and kinematics can improve our understanding of 
DLA systems (Wolfe et al. 1998; ); hydrodynamic simulations of merging protogalactic clumps are 



capable of explaining the large velocity spans seen in this class of absorption systems (Haehnelt 
et al. 1998), Here we have demonstrated the powerful tool which double LOSs can provide us. A 



convincing interpretation of our result, however, has to await a larger sample of such cases. 

Finally, let us emphasize that the scenario of small and compact C iv clouds surrounded by 
extended Ovi gas can indeed be theoretically understood and is in good agreement with the 



predictions of a cold-dark-matter based model (Rauch et al. 1997a ). Furthermore, in a hierarchical 
structure formation scenario damped Lya and Lyman limit systems can arise from relatively 
small (Mbaryon ~ 10 9 M ) merging protogalactic clumps. Even if the sizes derived using our 
photoionization models are not representative of transverse cloud sizes, the similar C IV line 
profiles and equivalent widths in A and B are still consistent with filamentary structures, which is 
also a prediction of hierarchical structure formation. Such a correlation between line profiles in A 
and B is, on the other hand, less evident for Mgn, implying gas inhomogeneities on spatial scales 
similar to the separation between LOSs, and in concordance with Mgn arising in smaller regions 
than C iv and O VI. 



S. L. thanks the BMBF (DARA) for support under grant No. 50OR96016; M. R. thanks 
NASA for support through grant HF-0107501-94A from the Space Telescope Science Institute; 
W. L. W. S. was supported by grant AST 9529073 from the National Science Foundation. A. S. 
thanks the financial support under grant no. 781-73-058 from the Netherlands Foundation for 
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Table 1. Journal of Observations. 



Object 


ID a 


UT date 


Coverage (A) 


S/N b FWHM (A) 


HE 1104-1805 A 


HST spectra 


Nov. 1995 


2222-3277 


20 


2.0 


HE 1104-1805 B 








17 




HE 1104-1805 A 


A AT spectra 


May 1993 


3170-7570 


82 


1.2 


HE 1104-1805 B 








30 




HE 1104-1805 A 


NTT spectra 


Feb. 1996 


3910-8290 


77 


0.8-1.0 


HE 1104-1805 B 








27 




HE 1104-1805 A 


Keck spectra 


Jan. and Feb. 1997 


3620-6080 


61 


0.08-0.13 


HE 1104-1805 B 








42 





& HST: Hubble Space Telescope; A AT: Anglo-Australian Telescope; NTT: New- Technology 
Telescope. 

b Maximum signal-to-noise ratio per pixel. 
c Reference: Paper I 
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Table 2. Lyman-edges in HE 1104-1805 A and B. 





^LLS 




log Ah i (cm 2 ) 


b (km s 1 ) 


A 


1.66162 




20.84 b 


27 




2.20016 


0.29±0.08 


16.67±0.12 


35 




2.29880 


0.14±0.07 


16.36±0.23 


40 


B 


1.66184 


2.34±0.51 


17.57±0.10 


30 




2.20099 


0.63±0.11 


17.00±0.08 


40 




2.29877 


0.21±0.11 


16.53±0.23 


35 



a Optical depth r = ln(/ + //_), where /+ is the 
extrapolated continuum redward of 912 x (1 + zlls) 5 an d 
/_ the absorbed continuum. 

b Column density constrained by the damped Lya wings. 
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Table 3. Line Parameters for z a b s — 1.66 Absorption Systems toward HE 1104-1805 A. 



Species 


Lines 


2 


Wi(A) 


7 /I — 1 \ 

b (km s J 


b 


log N (cm ) 


"logjV 


Spectr 


Hi 


1215 


1.661620 


19.05 


26.90 


0.30 


20.85 


0.01 


4 




1025 




2.99 










4 


















4 




918 














4 




1215 


1.664490 




28.50 


0.90 


16.59 


0.07 


4 




1025 














4 


















4 




918 














4 


Ci 


1656 


1.661620 


<0.02 


20.00 


0.00 


<12.80 




3 


Cn 


1036 


1.661640 




20.00 


0.00 


16.50 


0.20 


4 




1334 




0.46 










3 


Cm 


977 


1.661430 


1.62 


44.00 


0.00 


15.68 


0.20 


4 




977 


1.662800 




13.10 


0.00 


13.15 




4 




977 


1.664600 




15.80 


0.00 


15.42 




4 


Civ 


1548 


1.661186 


0.361 


15.88 


1.89 


13.31 


0.16 


1 




1550 




0.197 










1 




1548 


1.661412 




22.44 


3.14 


13.87 


0.06 


1 




1550 














1 




1548 


1.661694 




13.29 


2.94 


13.28 


0.21 


1 




1550 














1 




1548 


1.661908 




28.51 


5.09 


13.34 


0.13 


1 




1550 














1 




1548 


1.662776 


0.147 


15.11 


4.09 


13.13 


0.33 


1 




1550 




0.098 










1 




1548 


1.662856 




8.82 


0.69 


13.73 


0.07 


1 




1550 














1 




1548 


1.662984 




20.25 


12.40 


12.88 


0.55 


1 




1550 














1 




1548 


1.664327 


0.214 


4.69 


1.81 


12.32 


0.10 


1 




1550 




0.147 










1 




1548 


1.664598 




13.42 


0.94 


13.86 


0.03 


1 




1550 














1 




1548 


1.664739 




4.81 


0.97 


13.38 


0.08 


1 




1550 














1 




1548 


1.664860 




4.52 


0.68 


13.07 


0.03 


1 




1550 














1 


Ni 


1199 


1.661640 


0.29 


20.00 


0.00 


14.00 


0.20 


4 




1200.2 














4 




1200.7 














4 


Nil 


1083 


1.661640 


0.53 


20.00 


0.00 


15.00 


0.20 


4 


Nm 


989 


1.661430 


0.44 


44.00 


0.00 


13.60 


0.20 


4 


Nv 


1238 


1.662530 


<0.034 






<13.20 




4 


Oi 


988 


1.661640 




20.00 


0.00 


16.80 


0.20 


4 




1302 




0.75 










3 


Ovi 


1031 


1.662530 


0.93 


180.00 


0.00 


14.95 


0.20 


4 




1037 














4 


Mgi 


2852 


1.661523 


0.25 


14.20 


0.00 


12.33 


0.03 


2 




2852 


1.661793 




24.30 




11.65 


0.09 


2 


Mgn 


2796 


1.661370 


0.96 


3.61 


0.00 


10.44 


0.06 


2 




2803 




0.90 










2 




2796 


1.661549 




16.69 


0.00 


15.09 


0.03 


2 




2803 














2 
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Table 3 — Continued 



Species 


Lines 


z 


W r (A) a 


b (km s _1 ) 


<Tb h 


logiV(cm 2 ) 


<TlogiV 


Spcctr 




2796 


1.661710 




7.31 


0.00 


15.21 


0.18 


2 




zouo 














9 
Z 




97Qfi 






O. 1 u 


n nn 
u.uu 


1 1 ^n 


n n9 

u.uz 


2 




9&n^ 

ZOUO 














2 




97Qfi 




n^ 


i n nn 

-1U.UU 


n nn 
u.uu 


/ll on 
<v. 1 l.OU 




2 




9&n^ 

ZOUO 














2 




97Qfi 
z f yu 




u.uu 


7 7fl 


n nn 
u.uu 


12 19 


n n^ 
u.uo 


2 




zouo 














o 
Z 


Al I 


17^ 




<^0 09 
\U.UZ 


9fi nn 
zu.uu 




/io in 

<^ 1Z. 1U 






Alii 


1670 


1.661382 


0.404 


3.21 


0.47 


14.66 


0.31 


1 




1U I u 


I.UOIIOO 




12 29 


o.oo 


1 9 Q7 
iz. y ( 


n 1 7 

U. 1 l 






1U t u 


1 fifi1fiQ& 

l.uuiuyo 




1 9 77 

-1Z. ( I 


n Q9 

u.yz 


13 12 


n nfi 
u.uu 






1U I u 


i.uuiyyu 




7 1Q 

< . io 


n 


12 21 


n n9 

u.uz 


2 


Al in 


1 

lOOl 
1 &R9 


1 fifil ^78 
I.UOIO / o 


U. lOO 

n nQ9 
u.uyz 


9 QQ 
z.yy 


42 


12 34 


n n9 

u.uz 






1 S^zl 
lOO^t 


l.ODlOol 




1 Q 7Q 

iy . ( y 


1.01 


in oc; 

1Z.OO 


n n^ 

U.UO 


i 
i 




1 SR9 
lOOZ 














i 
i 




lOO^t 


1 fifil 797 
l.ODl f Z ( 




7 3A 


1 CiA 


1 9 97 
IZ. Z i 


n no 
u.uy 


i 
i 




1 &R9 


















1854 


1.661981 




5.32 


1.93 


11.58 


0.08 


1 




1862 














1 


Si i 


1845 


1.661620 


<0.01 


20.00 




<12.90 




1 


9l TT 
Ol 11 


i &n& 

louo 
1 ^9^ 

lOZU 


1 fifil ^71 
l.tJlJlO I 1 


U. -1UO 

n A^n 


fii 

O.U1 


n i q 
u.iy 


1 ^ 1 

IO. io 


n n9 

u.uz 


2 
j 




i &n& 

louo 


l.UUlOOO 




1 fi fiQ 

-iu.uy 


n 

u.oy 


1 A Q7 


n m 

U.Ul 






1526 














1 




louo 


1 fifil 77^ 
l.UOl 1 I o 




7 TI 

I .Ol 


1 n& 

l.UO 


1 1 Qfi 

io. yu 


n no 
u.uy 


2 




1526 














1 




1808 


1.661984 




8.70 


0.59 


13.66 


0.02 


1 




lOZU 














2 


9i TTT 
Ol 111 


1ZUU 


l.UU14:OU 


n fin 
u.uu 


44 00 


n nn 
u.uu 


1 a nn 

I'i.UU 








1ZUU 


l.UU'iUUU 




1 ^ £n 


n nn 
u.uu 


1 ^ 7n 

IO. < u 


n 9n 




Ol IV 


1 ^Q^ 

loyo 
1402 
loyo 
1402 


l.UUlU'iU 

1 Rfi9snn 


fiQ 
\u.oy 

^n ^a 
m 7 






<<1 a ns 




2 

2 

i 
i 

2 




loyo 


l.OO^tOOU 


n 1 7n 

U. 1 I u 






10.11 


u.u* 


i 
i 




1402 




0.044 










i 


Ti II 


1910.6 
1910.9 


l.DDlooy 


0.007 
0.003 


2.81 


1.16 


12.28 


0.05 


i 
i 


Cr II 


2056 
2062 
2066 


1.661381 


0.051 
0.043 
0.028 


3.94 


0.22 


12.86 


0.01 


i 
i 
i 




2056 


1.661575 




8.26 


0.73 


12.61 


0.02 


i 




2062 














i 




2066 














i 




2056 


1.661705 




1.75 


0.74 


11.68 


0.12 


i 




2062 














i 




2066 














i 


Mm 


2795 


1.661420 


0.06 


14.20 


0.00 


11.80 


0.20 


2 


Mnn 


2576 
2594 


1.661469 


0.05 


14.20 


0.00 


12.66 


0.06 


2 
2 
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Table 3 — Continued 



Species 


Lines 


z 


VK r (A) a 


6 (km s- 1 ) 


<r b b 


log-/V(cm~ 2 ) 


"logjV 


Spectrum 




2606 














2 




2576 


1.661734 




24.30 




11.82 


0.38 


2 




2594 














2 




2606 














2 


Fe I 


2484 


1.661620 


< 0.02 


20.00 




<12.30 




2 


Fe II 


1608 


1 ^/"l n "7/" 

1.661376 


0.280 


4.80 


0.16 


14.49 


0.01 


1 




2249 




0.045 










1 




2260 




0.079 










1 




1608 


1.661549 




8.49 


0.70 


14.29 


0.03 


1 




2249 














1 




2260 














1 




1608 


1.661709 




8.83 


1.11 


13.85 


0.06 


1 




2249 














1 




2260 














1 




1608 


1.661983 




9.12 


1.10 


13.34 


0.04 


1 




2249 














1 




2260 














1 


Nin 


1709 


1.661382 


0.052 


3.05 


0.35 


13.07 


0.02 


1 




1741 




0.051 










1 




1751 




0.048 










1 




1709 


1.661559 




11.52 


1.38 


12.98 


0.04 


1 




1741 














1 




1751 














1 




1709 


1.661716 




0.52 


0.19 


12.29 


0.27 


1 




1741 














1 




1751 














1 




1709 


1.661921 




21.67 


7.71 


12.47 


0.12 


1 




1741 














1 




1751 














1 


Znll 


2062 


1.661373 


0.026 


3.04 


0.24 


12.37 


0.01 


1 




2026 




0.042 










1 




2062 


1.661560 




7.64 


1.21 


11.82 


0.04 


1 




2026 














1 



a Rest-frame equivalent widths. For blends, the total equivalent width is given. Upper limits (3<r) 
represent non-detections. 

b <T[, = indicates fixed b parameter. 

c (l) Keck; (2) NTT; (3) AAT; (4) HST 
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Table 4. Ionic Column Densities (cm 2 ) and Gas-Phase Abundances in the z a b s 

toward HE 1104-1805 A. 



1.66162 DLA 



Species 


logJV flt (Xi) 




logJV a pp(X<) 


"logiV 


log7V ad (Xi) 


°"IogiV 


log TV (X) 




[X/Hf 


ff[X/H] 


U T 

11 1 


ZU.oD 


n ni 

U.Ul 






ZU.oD 


U.Ul 










P T 
^ 1 


/I') Qfl 
s,lZ.OU 








/io on 




1 « ^0 
1D.OU 


ft 9ft 


n qi 
-u.yi 


n 9n 
u.zu 


Cn 


16.50 


0.20 






16.50 


0.20 










Cm 


15.70 


0.20 






15.70 


0.20 










Civ 


14.13 


0.05 


14.13 


0.01 


14.13 


0.05 










N I 


14.00 


0.20 






14.00 


0.20 


15.04 


0.16 


-1.86 


0.16 


Nil 


15.00 


0.20 






15.00 


0.20 










Nm 


13.60 


0.20 






13.60 


0.20 










Nv 


<13.20 








<13.20 












Oi 


16.80 


0.20 






16.80 


0.20 


16.80 


0.20 


-0.98 


0.20 


Ovi 


14.95 


0.20 






14.95 


0.20 










Mgi 


12.41 


0.09 






12.41 


0.09 


>15.46 


0.09 


>-0.97 


0.09 


Mgn 


15.46 


0.09 






>15.46 


0.09 










All 


<12.10 








<12.10 




14.56 


0.15 


-0.77 


0.15 


Alii 


14.55 


0.24 


13.43 


0.16 


14.55 


0.24 










Aim 


13.06 


0.02 


13.08 


0.01 


13.06 


0.02 










Si i 


<12.90 








<12.90 




15.38 


0.02 


-1.02 


0.02 


Sin 


15.38 


0.02 


15.26 


0.01 


15.38 


0.02 










Sim 


14.00 


0.20 






14.00 


0.20 










Si iv 






<14.08 




<14.08 












Tin 


12.28 


0.05 


12.28 


0.07 


12.28 


0.07 


12.28 


0.07 


-1.50 


0.07 


Cm 


13.07 


0.01 


13.05 


0.01 


13.07 


0.01 


13.07 


0.01 


-1.46 


0.01 


Mm 


11.80 


0.20 






11.80 


0.20 


12.77 


0.14 


-1.61 


0.14 


Mnii 


12.72 


0.18 






12.72 


0.18 










Fei 


<12.30 








<12.30 




14.77 


0.01 


-1.59 


0.02 


Fen 


14.77 


0.01 


14.68 


0.02 


14.77 


0.02 










Nin 


13.42 


0.03 


13.43 


0.02 


13.42 


0.03 


13.42 


0.03 


-1.68 


0.03 


Znll 


12.48 


0.01 


12.45 


0.01 


12.48 


0.01 


12.48 


0.01 


-1.02 


0.01 



^Relative to solar element abundances compiled by Verner et al. 1994 
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Table 5. Line Parameters for the z a ^ s — 1.66 Absorption Systems toward HE 1104-1805 B. 



Species 


Lines 


z 


T T 7 ( % \ a 
W r (A) 


b (km s ) 


b 


1 AT / — 2 \ 

log N (cm ) 


flogJV 


Spcctr 


Hi 


1215 


1.661840 


1.690 


30.00 


0.00 


17.57 


0.00 


4 




1025 




0.790 










4 




918 














4 
4 




1215 


1.664970 




30.00 


0.00 


15.90 


0.13 


4 




1025 














4 




918 














4 
4 


Ci 


1656 


1.661650 


<0.09 


20.00 


0.00 


<13.50 




3 


Cn 


1036 


1.661650 




20.00 


0.00 


14.50 


0.20 


4 




1334 




0.18 










3 


Cm 


977 


1.664480 


0.95 


37.00 


0.00 


15.50 


0.20 


4 




977 


1.661830 




29.30 


0.00 


14.00 


0.20 


4 


Civ 


1548 


1.661130 


0.630 


19.93 


2.33 


13.24 


0.04 


1 




1550 




0.361 










1 




1548 


1.661512 




8.82 


0.88 


13.56 


0.06 


1 




1550 














1 




1548 


1.661759 




21.96 


4.49 


13.81 


0.10 


1 




1550 














1 




1548 


1.662036 




12.62 


1.74 


13.77 


0.09 


1 




1550 














1 




1548 


1.662250 




47.39 


4.03 


13.71 


0.15 


1 




1550 














1 




1548 


1.662878 




20.75 


3.44 


13.31 


0.10 


1 




1550 














1 




1548 


1.663167 




5.89 


2.53 


12.38 


0.16 


1 




1550 














1 




1548 


1.663384 




8.02 


2.11 


12.55 


0.08 


1 




1550 














1 




1548 


1.664721 


0.077 


3.54 


1.49 


12.61 


0.07 


1 




1550 




0.067 










1 




1548 


1.664887 




5.43 


0.36 


13.63 


0.03 


1 




1550 














1 


Ni 


1199 


1.661650 


<0.11 


20.00 




<14.00 




4 


Nil 


1083 


1.661650 


<0.13 


20.00 




<13.70 




4 


Nm 


989 


1.661830 


0.27 


37.00 


0.00 


14.00 


0.20 


4 




989 


1.662840 




29.90 


0.00 


13.70 


0.20 


4 




989 


1.665070 


0.27 


37.20 


0.00 


14.20 


0.20 


4 


Nv 


1238 


1.662530 


<0.19 






<14.00 




4 


Oi 


988 


1.661650 


<0.13 


20.00 




<14.10 




4 




1302 




<0.19 










3 


Ovi 


1031 


1.662530 


0.86 


110.00 


0.00 


14.95 


0.20 


4 




1037 














4 


Mgi 


2853 


1.661650 


<0.10 


20.00 




<11.90 




2 


Mgn 


2796 


1.661557 


0.343 


7.83 


0.00 


13.14 


0.13 


2 




2803 




0.239 










2 




2796 


1.661783 




26.52 


0.00 


12.73 


0.09 


2 




2803 














2 




2796 


1.664930 


<0.09 


10.00 




<12.40 




2 




2803 














2 


All 


1765 


1.661650 


<0.10 


20.00 




<12.60 




1 


Alii 


1670 


1.661556 


0.097 


7.70 


0.00 


12.26 


0.03 


1 
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Table 5 — Continued 



Species 


Lines 


2 


Wr (A) a 


b (km s 1 ) 




log iV(cm~ 2 ) 


ClogiV 


Spectrum 




1670 


1.661781 




14.10 


0.00 


11.80 


0.06 


1 




1670 


1.661990 




8.10 


0.00 


11.45 


0.10 


1 


Aim 


1854 


1.661564 


0.041 


7.75 


1.29 


12.08 


0.11 


1 




1862 




0.036 










1 




1854 


1.661784 




14.07 


13.05 


11.98 


0.35 


1 




1862 














1 




1854 


1.661956 




8.15 


11.22 


11.81 


0.46 


1 




1862 














1 


Si I 


1845 


1.661650 


<0.08 


20.00 




<13.20 




1 


Sin 


1526 


1.661552 


0.076 


8.30 


0.94 


13.29 


0.03 


1 




1526 


1.661747 




7.72 


3.94 


12.80 


0.15 


1 




1526 


1.661941 




8.99 


3.69 


12.89 


0.12 


1 


Sim 


1206 


1.661830 


0.55 


37.00 




14.10 


0.20 


4 


Si iv 


1393 


1.661840 


1.166 






<14.28 




1 




1402 




0.580 










1 




1393 


1.664930 








<12.86 




1 




1402 




<0.033 










1 


Fen 


2344 


1.661578 


0.012 


7.86 


4.30 


12.47 


0.08 


2 




2382 




0.030 










2 




2600 




0.045 










2 



a Rest-frame equivalent widths. For blends, the total equivalent width is given. Upper limits (3<r) 
represent non-detections. 

b cr f) = indicates fixed b parameter. 

c (l) Keck; (2) NTT; (3) A AT; (4) HST 
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Table 6. Ionic Column Densities (cm~ 2 ) in the LLS at z ahs = 1.66184 toward HE 1104-1805 B in 
Comparison to Results from Photoionization Models. 



Species 


logiVfitpQ) 


flogJV 


logJVappCXi) 


CT logJV 


log AT ad pQ) 


""logJV 


MODEL l a 


MODEL 2 b 


Hi 


17.57 


0.10 






17.57 


0.10 


17.57 


17.57 


Hii 














19.89 


19.82 


Ci 


<13.50 








<13.50 




11.87 


12.14 


Cn 


14.50 


0.20 






14.50 


0.20 


14.47 


14.75 


Cm 


15.50 


0.20 






15.50 


0.20 


15.25 


15.24 


Civ 


14.41 


0.05 


14.40 


0.01 


14.00 c 


0.06 


13.99 


13.99 


Ni 


<14.00 








<14.00 




10.55 


10.97 


Nil 


< 13.70 








< 13.70 




13.00 


13.31 


Nm 


14.00 


0.20 






14.00 


0.20 


13.79 


13.77 


Nv 


<14.00 








<14.00 




11.30 


10.87 


Oi 


<14.10 








<14.10 




12.52 


12.84 


Ovi 


14.95 


0.20 






14.95 


0.20 


11.54 


10.76 


Mgi 


<11.90 








<11.90 




11.61 


11.68 


Mgn 


13.28 


0.09 






13.28 


0.09 


13.23 


13.29 


All 


<12.60 








<12.60 




10.69 


11.00 


Alii 


12.44 


0.03 


12.44 


0.01 


12.44 


0.03 


12.98 


13.26 


Aim 


12.45 


0.05 


12.39 


0.03 


12.45 


0.05 


12.60 


12.86 


Si i 


<13.20 








<13.20 




10.29 


10.48 


Sin 


13.52 


0.05 


13.52 


0.01 


13.52 


0.05 


13.30 


13.47 


Sim 


14.10 


0.20 






14.10 


0.20 


14.06 


14.12 


Si iv 






<14.28 




<14.28 




13.43 


13.25 


Fen 


12.47 


0.08 


11.78 


0.12 


12.47 


0.08 


11.45 


11.08 



a MODEL 1: Haardt & Madau, logT = -2.95, logn H = -1.8, Z = 0.63 Z DLA , S = 1.6 kpc. 
b MODEL 2: Power Law a = -2.0, logT = -3.02, logra H = -1.91, Z = 0.79 Z DLA , S = 1.8 kpc. 
c Sum of fit components 2 and 3. 
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Table 7. Ionic Column Densities (cm 2 ) in the Civ Systems at z a b s = 1.66280 toward HE 

1104-1805 A and Photoionization Models. 



Species 


log N Rt {Xi) 


°"logJV 


logiVapppQ) 


"logiV 


logAT ad pQ) 


"■logJV 


MODEL a 


Hi 


16.59 


0.07 






16.05 b 


0.07 


16.05 


Hii 














18.92 


Cm 


13.15 


0.20 






13.15 


0.20 


14.42 


Civ 


13.87 


0.10 


13.85 


0.01 


13.87 


0.10 


13.84 


Niii 


<13.00 








<13.00 




12.97 


Mgn 


<11.80 








<11.80 




11.76 


Si iv 






<12.58 




<12.58 




12.88 



a MODEL: Haardt & Madau, logT = -2.45, logn H = -2.30, Z = Z DLA , S = 0.5 kpc. 
b A r (Hl) distributed according to the TV(Mgll) ratios. 
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Table 8. Ionic Column Densities (cm 2 ) in the Civ Systems at z a b s — 1.66465 toward HE 

1104-1805 A and Photoionization Models. 



Species 


logATfitpQ) 


°"logjV 


logiVapppQ) 


°"logiV 


logAT ad pQ) 


ClogJV 


MODEL a 


Hi 


16.59 


0.07 






16.44 b 


0.07 


16.44 


Hii 














19.22 


Cm 


15.42 


0.20 






15.42 


0.20 


14.74 


Civ 


14.04 


0.03 


14.03 


0.01 


14.04 


0.03 


14.06 


Nm 


<13.20 








<13.20 




13.29 


Mgn 


12.19 


0.03 






12.19 


0.03 


12.20 


Sim 


13.70 


0.20 






13.70 


0.20 


13.42 


Si iv 






13.11 


0.07 


<13.11 


0.07 


13.19 



a MODEL: Haardt & Madau, logT = -2.54, logn H = -2.21, Z = Z DLA , S = 0.9 kpc. 
b W(Hl) distributed according to the TV(Mgll) ratios. 
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Table 9. Ionic Column Densities (cm 2 ) in the Civ Systems at z a b s = 1.66493 toward HE 

1104-1805 B and Photoionization Models. 



Species log iV flt (X^ ) a logN log iY ap p (Xi ) <Ti og jv log iV ad pQ) <r logJV MODEL a 



Hi 


15.90 




15.90 


15.90 


Hn 








18.90 


Cm 


<12.90 




<12.90 


14.17 


Civ 


13.67 0.03 


13.61 


0.02 13.67 0.03 


13.67 


Niii 


<13.20 




<13.20 


12.71 


Mg II 


<12.40 




<12.40 


11.33 


Si IV 




<12.86 


<12.86 


12.58 


a MODEL: Haardt & Madau, logT = 


—2.35, logriH 


= -2.40, Z = 0.63 Z DLA , S = 


0.6 kpc. 
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Fig. 1. — HST FOS spectra of HE 1104-1805 A and B and their la errors. Fluxes have been 
corrected for galactic extinction (Seaton |1979| ) with E(B — V)=0.09 (Reimers et al. 1995). Also 



shown are the AAT spectra which have been scaled to the HST flux levels, smoothed and rebinned 
to the FOS resolution, and plotted for A > 3277 A. Metal absorption lines lying longward of Lya 
emission were covered by the NTT and Keck spectra. The power law and the QSO continuum 
derived as described in the text have been overplotted. Expected positions of emission lines by 
Lya, Ly/3, Nv, Ovi, and Svi are indicated. 

Fig. 2.— Sections of the normalized HST FOS (FWHM ~ 2 A) and AAT (FWHM ~ 1.2 A) 
spectra of HE 1104-1805 A (upper panels) and B, showing synthetic Voigt profiles of the high- 
order Hi Lyman-series lines (top) and most relevant metal lines belonging to the damped Lya 
system (A) and Lyman- limit System (B) at z = 1.66. Tickmarks (from upper panel to right-lower 
panel) indicate the positions of Hi AA918 to 949, Cn AA1036, 1334, Nil A1083 and Oi AA988, 1302. 



Fig. 3. — Low ions: Velocity profiles at 6.6 km s resolution of lines associated with the damped 
Lya system toward HE 1104-1805 A (left-hand panel), and the Lyman-limit System toward HE 
1104-1805 B. The Mgn line profiles are at 40 km s _1 resolution. The zero velocity scale in both 
panels corresponds to z = 1.66164. The smoothed lines represent the best-fit Voigt profiles and 
ticks mark the line positions. 

Fig. 4. — High ions: Velocity profiles at 6.6 km s _1 resolution of lines associated with the damped 
Lya system toward HE 1104-1805 A (left-hand panel), and the Lyman-limit System toward HE 
1104-1805 B. The Al II A1670 velocity profiles are also shown for comparison purposes (see text). 
The zero velocity scale in both panels corresponds to z = 1.66164. The smoothed lines represent 
the best-fit Voigt profiles and ticks mark the line positions. The absorption feature at v ~ +110 
km s" 1 in the Civ A1548 B plot is identified with Sin A1808 at z = 1.28. 

Fig. 5. — Left: Velocity profiles at 6.6 km s _1 resolution of Cm, Znii and Tin lines associated 
with the damped Lya system toward of HE 1104-1805 A. Right: Expected positions of the same 
ions in HE 1104-1805 B. The zero velocity scale in both panels corresponds to z = 1.66164. The 
smoothed lines represent the best-fit Voigt profiles and ticks mark the line positions. 

Fig. 6. — Bottom: Velocity profiles of the Civ AA1548,1550 (at 6.6 km s _1 resolution) and Mgn 
AA2796,2803 (at 40 km s _1 resolution) doublets in A (left) and B. In both panels v = km s _1 
corresponds to z = 1.66164. In the Mgn A2796 plot, the feature at —100 km s _1 is identified with 
Mm A2795. Top: HST velocity profiles of the corresponding Lya, Ly/3 and Ly7 Hi lines (note the 
different velocity scales). 

Fig. 7. — Section of the HST spectra A and B showing the O vi AA1031,1037 doublet at z = 
1.66253. Solid lines: Voigt profiles — convolved with a FWHM = 2 A Gaussian — with b = 180 
km s _1 (A) and 110 km s _1 (B), and column density log N = 14.95. Dotted lines: Voigt profiles of 
3 (A) and 2 (B) O vi lines with common b = 50 km s _1 and total column density log N = 15.1 (A) 
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and 15.0 (B). The lines are placed at redshifts labeled 1, 2 and 3 (A), and 4 and 5 (B) in Fig. ^. 

Fig. 8. — Section of the AAT spectra A and B. Tickmarks indicate the positions at which the N v 
AA1238, 1242 doublet lines would appear if N v were present at z = 1.66253. Part of the red wing 
of the damped Lya line appears in the overplotted synthetic spectrum A. In B, the absorption at 
A = 3308 A is identified with a Lya line, supported by the detection of Ly/3 at the same redshift 
in the HST spectrum. 

Fig. 9. — Hydrogen column density vs. metallicity in the O vi-phase for different CLOUDY models 
assuming the Haardt & Madau continuum (solid line) and a power law as ionizing background. 
The relative abundances are the same as in the DLA gas. These models yield throughout 
logJV(Ovi)= 15.0 and logJV(Nv)= 13.2. 

Fig. 10. — Column densities vs. ionization parameter T in the O vi-phase for different CLOUDY 
models. Solid line: Haardt & Madau radiation field as ionizing background, and Z = Z-qla- Dotted 
line: Same radiation field, with Z = 0.4 Zdla- Dashed line: Power law and Z = Zdla- The relative 
abundances are the same as in the DLA gas. These models yield throughout logiV(Ovi)= 15.0. 
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